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The composition-independent virial coefficients of a d-dimensional binary mixture of (additive) 
hard hyperspheres following from a recent proposal for the equation of state of the mixture [Santos, 
A., Yuste, S. B., and Lopez de Haro, M., f999, Molec. Phys., 96, 1] are examined. Good 
agreement between theoretical estimates and available exact or numerical results is found for d — 2, 
3, 4 and 5, except for mixtures whose components are very disparate in size. A slight modification 
that remedies this deficiency is introduced and the resummation of the associated virial series is 
carried out, leading to a new proposal for the equation of state. The case of binary hard-sphere 
mixtures (d = 3) is analyzed in some detail. 



Depending on the nature of the independent variables, it is well known that the full thermodynamic description 
of a given system requires the availability of either one of the thermodynamic potentials or of the corresponding 
equations of state. In the case of classical fluids, the usual expression for the equation of state (EOS) is a relationship 
among the pressure p, the density p and the temperature T of the fluid. However, except for very few ideal systems, 
the explicit (exact) form of this relationship is not known in general. Therefore, given the fact that an experimental 
determination of the EOS for every particular fluid is of course not practical, researchers have concentrated their 
efforts in proposing "reasonable" empirical or semi-empirical approximations. Perhaps the neatest example of such 
an approach, already produced over a century ago and providing an essentially correct qualitative picture of the 
thermodynamic properties of gases and liquids, is the celebrated van der Waals EOS An alternative proposal, first 
introduced by Kammerlingh Onnes j^j as a mathematical representation of experimental results on the equation of 
state of gases and liquids, is the so called virial equation of state. Written for the compressibility factor Z = pj pksT 
(where ks is the Boltzmann constant) as a power series in p, it provides an expansion that gives the deviation from 
ideal gas behaviour in ascending powers of the density, namely 



Here, the B n are the virial coefficients. For a simple fluid, the virial coefficients are in general functions of T alone, 
while in the case of mixtures they also depend on composition ||. An interesting aspect of the virial equation of 
state, discovered many years after it was introduced, is the fact that it can be derived rigorously using statistical 
mechanics This in turn implies that the virial coefficients are not merely empirical constants but are rather 
related to intcrmolccular interactions in a well-defined manner. It is unfortunate, however, that in general the actual 
computation of the virial coefficients is a formidable task and that the radius of convergence of the series in Eq. (|l|) 
is not known. 

In the case of hard-core fluids (rods, discs, spheres and hyperspheres) the virial coefficients are easier to compute than 
those corresponding to any other intermolecular interactions. Even in these model systems and with the exception of 
hard rods, where the exact EOS and corresponding virial coefficients are known, only the first few of such coefficients 
are available so far ||-|| because the number of cluster integrals involved increases very rapidly with the order of the 
coefficient. The availability is scarcer for mixtures than for simple fluids, although recently there has been renewed 
interest in partially remedying this deficiency [p|-^9[ . The limited knowledge of the virial coefficients has proved useful, 
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however, in the sense that proposals for the EOS of these systems may be judged amongst other things by how well 
the virial coefficients arising in such proposals compare with the exact values. The rationale here is that a "good" 
theoretical EOS should lead to an accurate prediction of the value of the virial coefficients. Alternatively the few 
available virial coefficients may also be profitably used to construct a rational (Pade or Levin) |]^,po|-p^| aproximation 
to the EOS or to guide in the construction of a theoretical EOS, as was the case with the Carnahan-Starling (CS) 
EOS for a simple hard-sphere fluid. 

A varietyof (approximate) EOS for simple fluids composed of hard discs and hard spheres are available in the 
literature and some work has also been reported on simple fluids of hard hyperspheres [p7|-|3~i| . Very recently we 
introduced a simple recipe to derive the compressibility factor of a multicomponent mixture of c?-dimensional hard 
hyperspheres in terms of that of the corresponding single component system [32). A straightforward consequence of 
such a recipe is that one can easily derive the explicit expressions for the virial coefficients of the mixture in terms 
of those of the single component fluid for all d. The major aim of this paper is to further assess the usefulness of 
our recipe by comparing the predictions of the values of the virial coefficients of binary mixtures of hard discs, hard 
spheres and hard hyperspheres in d = 4 and d = 5 with those obtained through Monte Carlo integration. Along 
the way we will introduce slight modifications to our original proposal and consider an alternative EOS to cope with 
mixtures very disparate in size. This new EOS is consistent with the forms suggested by Wheatley for the virial 
coefficients of hard discs |53| and hard spheres |34| . 

The paper is organized as follows. In the next Section and in order to make the paper self-contained, we will briefly 
recall the main results of our previous work and write down explicitly the composition-independent virial coefficients 
of a d-dimensional binary mixture of (additive) hard hyperspheres. Section III deals with a new proposal for the EOS 
of the binary mixture that attempts to incorporate additional features not present in the original recipe. In Section 
IV we focus on the case of hard spheres because of its intrinsic interest. We close the paper in Section [v| with an 
assessment of the different results. Some mathematical details of our derivations are presented in two Appendices. 



II. APPROXIMATE EQUATION OF STATE AND VIRIAL COEFFICIENTS FOR A BINARY MIXTURE 

OF d-DIMENSIONAL HARD SPHERES 



In this Section we first consider a multicomponent mixture of d-dimensional (additive) hard hyperspheres and then 
restrict ourselves to binary mixtures of hard hyperspheres in d dimensions. 

Let the number of components be N, the total number density of the mixture be p, the set of mole fractions be 
{xi}, and the set of diameters be {<Ji}, (i — 1,2, ... , N). The packing fraction is r\ — rji = Vdp(c d ), where r\i — 

VdPi&f is the partial packing fraction due to species i, pi — pxi is the partial number density corresponding to species 
i, Vd = (7r/4) d / 2 /r(l + d/2) is the volume of a d-dimensional sphere of unit diameter and (<j™) = YliLi x i (J i- I n 
previous work |32| we proposed a simple equation of state for the mixture, Z m (ry), consistent with a given EOS for a 
single component system, Z s (j]), at the same packing fraction rj, namely 



Z m (rj) = 1 + [Z s (r,) - 1] 2 1 " d A 



(2) 



with 



i^+p- 1 ) g {d+ P -iy. 



rd\2 



m\(d + p — 1 — m)! 



(3) 



In the one-dimensional case, Eq. (g) yields the exact result Z m (i]) = Z s (rj). Further, for binary mixtures in d — 2, 3, 
4 and 5, it proved to be very satisfactory when a reasonably accurate Z s (rj) was taken [ p4] , ^2|j3^j3^ | . 

From the virial expansion Z s (rj) = 1 + X)^°=2 k™ 7 ?™ ' wnere b n are (reduced) virial coefficients, and according to 
Eqs. (0)-(0), the n-th virial coefficient of the mixture is given by 
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and explicitly for n > 2 
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Now we particularize to binary mixtures. In that case, it is useful to define composition-independent coefficients 
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Our objective is to use our model to get approximate expressions of B nin2 as explicit functions of <7i, er 2 , and d. In 
Appendix |X| it is shown that the result is 
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The special case of n = 2 is obtained by first setting b n = 2 d ~ 1 and then setting n = 2. The result is 

C&U = 2 d - 2 n 1 (ni - 1), ^=0, CW ina = inina. (11) 

This gives the exact second virial coefficient. 
Equation ([?]) expresses the reduced coefficient 

B* un2 (a) = B nun2 {a l ,a 2 )a- d{n - 1) a- d ^ (12) 

as a polynomial in a of degree d. Although it is known that the exact coefficients do not have in general such a 
polynomial structure |$4j , most of the proposals to date for d = 2 and d = 3 are polynomials ]T^ , |l6| , ^8| , p3| , p4||3^j3^] . 
Let us see which consistency conditions the approximation (0)-(^) fulfills. First, it verifies the obvious symmetry 
property B ni>n2 {ai, er 2 ) = B n . 2ni (<r 2 , <7i), which implies that B^ i n2 (a) = a d B^ 2 ni (l/a). In addition, the coefficients 
(||)-(|of) satisfy the properties 
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Equation (Il3|) implies that 
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i.e., if both species have the same size we recover the one-component case. The same situation appears if x± = or 
x 2 = 0, which means that 23* (a) = v'^~ 1 b n a d and 23g „(a) = v , j~ 1 b n . This is verified as a consequence of Eq. (|l5|). 
A subtler consistency condition is pifl 
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which ensures that the derivative of the excess free energy of mixing with respect to a 2 is zero when the spheres are 
of the same size [M . In the approximation (j?]) , Eq. ( |l8| ) requires 



dCi% 2 +(d- l)Cr« na + 2 d -\d - l)Ci% 2 + 2 d -\d + l)C&\ ni + Cg 2 ] ni = <2^6 n 



(19) 

If we interchange the roles of n\ and n.2 in Eq. ( |19[) and add the original and the transformed equations, we get Eq. 
(pf). If, instead, we subtract both equations, we obtain 
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Thus, enforcement of Eq. (|T^) is equivalent to enforcement of Eq. (p0f). The property ( |i"4| ) guarantees that the 
condition @ is fulfilled by the coefficients (I)-©. 

The most stringent conditions appear in the limits a — ► and a — > 00. In that case, Eq. (Q) yields 
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On the other hand, the exact result is [B3|34l3 
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The above condition means that in the limit of infinite size asymmetry, the smaller spheres contribute to the total 
pressure as if they were a one-component system in a free volume equal to the total volume minus the volume occupied 
by the larger spheres. Our approximation (0)-([n]) gives the correct forms for the asymptotic behaviours of 23* „ 2 , 
but not the correct coefficients (except for n 2 — 1 in the limit a — > and for m = 1 in the limit a — > 00, apart from the 
trivial cases n\ — 0, n and n 2 — 0, n). In fact, in our approximation the virial coefficient 23„ depends parametrically 
on b n only, and not on the previous coefficients b ni ,n\ < n, as is needed in Eq. (|2^). This is a consequence of the form 
of our approximation in which the compressibility factor of the mixture at a given packing fraction 77 is expressed in 
terms of the compressibility factor of the monodisperse system at the same packing fraction. 
Another consistency condition is flij 
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This comes from the exact value of the derivative of the excess free energy with respect to molecular size in the limit 
of infinite size asymmetry |38|. In general, this condition is not verified by our approximation either. We consider 
the condition (|2|) to be more important than ( |23| ) because it is not restricted to n% = 1 (or n\ = 1) and also because 
it refers to the leading term in the limit of infinite size asymmetry, while the condition (p3) refers to the sub-leading 
term. 



TABLE I. Reduced virial coefficients b n for the one-component case. 
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b Refs. 



3.12801775 
4.25785446 
5.336897 
6.3626 
7.351 
8.338 



4 
10 

18.36477 
28.2245 
39.739 
53.539 
70.78 



32.4057594 
77.7451797 

145.9 

252.0 
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311.18341 
843.4 
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Table | gives the known values of the coefficients b n for 2 < d < 5. From them, Eqs. (^-(p"l|) can be used to 
estimate the coefficients £?* for any value of a in these dimensions. We will come back to this point later on. The 
values of Table | can also be used to check to which degree Eq. ( p2| ) is violated by the approximation (0)-(^O|). The 
less favorable case corresponds to n\ = 1 in the limit a — > (or, equivalently, n% = 1 in the limit a — > oo). Let us 
then define 



rt(3) , r (i) 
(n - l)/n ]b n -i 



According to Eqs. (||) and we have 



_ (2-' + .-3)l + 2<--l 
(n - l)6„_i 

Table || shows the values of R n that can be obtained from the known virial coefficients. In the cases d = 2 and 
d = 3 the known values of the ratio R n differ from 1 less than 5%. In fact, by assuming R n ~ 1, Eq. (|2^) can be 
used to estimate 6 9 and 6 10 . This gives b 9 ~ 9.39 (d = 2), 6 9 ~ 86.7 (d = 3), 6 10 ~ 10.44 (d = 2) and b w ~ 103.6 
(d = 3). Estimates based on Pade approximants are [|| 69 ~ 9.37 (d = 2), 69 ~ 93.1 (d = 3), 610 — 10.55 (d = 2) and 
610 — 123.2 (d = 3). On the other hand, the deviations of R n from 1 become more important as the dimensionality 
increases and are already relatively large for d = 5. Thus, at least for those high dimensionalities, it appears to be 
convenient to modify the prescription (Q) as to satisfy the requirement (|22]). This is the subject of the next Section, 
where we will also discuss in more detail how good the estimates of B* t , as obtained after using Eqs. (|)-® with 
the values of the coefficients b n given in Table |, are for any value of a and 2 < d < 5. 

TABLE II. Ratio R n , as given by Eq. @. 

n d = 2 d= 3 d = 4 d = 5 

3 1.032 1 0.944 0.883 

4 1.014 1.018 1.072 1.148 

5 0.999 1.001 1.078 1.452 

6 0.991 1.007 1.132 0.736 

7 0.989 1.023 

8 0.992 1.047 
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III. A NEW PROPOSAL FOR THE EQUATION OF STATE OF A BINARY MIXTURE 



In order to account for the conditions imposed by Eq. fl22| ) without greatly sacrificing the simplicity of our original 
recipe, let us now explore the possibility of keeping the structure (Q), but with expressions for the coefficients different 
from The requirement ( p2| ) implies that 
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In addition, the condition Jl7j ) yields Eq. Jl6|). Note that Eq. ( [l3| ) is more restrictive than Eq. (|l6|). Since we want 
our modified expressions to remain as close as possible to the original ones, we will take @ as a second condition. 
From Eqs. (|26|) and (Of) we have 
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We need a third condition to close the problem. Since Cn},n 2 does not enter into Eq. (|l]) and in agreement with our 
philosophy of departing from the original formulation as less as possible, one possibility is to keep Eq. (^); however, 
the resulting approximation is consistent with the conditions ( [l5| ) but does not satisfy the requirement ([l9|). If we 
enforce the fulfillment of Eq. ( |l9| ) [or, equivalently, of Eq. (p(j|)], then Cni,n 2 must have the following form 
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where -D nii „ 2 must be symmetric, but otherwise remains so far unknown. To close the problem, we impose the more 
restrictive property (|l4|). This gives rise to 
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Equations (|27D-(|30D close the modified version of our approximation. 

It must be noted that Eqs. (B9 ) and (pOl) are meaningless if d = 2. This is related to the fact that if d = 2 it is 
impossible to enforce conditions (jl7|), (|lq) and ( |22] ) simultaneously on a polynomial of second degree for B* n , as 



the number of conditions exceeds the number of unknowns. As a matter of fact, Eqs. (27) and (E8[) are enough to 
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uniquely determine £?* „ 2 in the case d — 2 7 irrespective of the value of Cni,n 2 - The result is 
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This was the form proposed by Wheatley ]33| by imposing conditions ([l7]) and (|2^). However, condition (|T^) is not 
verified unless ni(b n — b ni ) = n 2 {b n — b n2 ), which is only true if b n — n, i.e., if one considers the scaled particle theory 
p9[ . For the sake of comparison, the form that stems from our original approximation, Eqs. (f7])~(fIo|), is (for d = 2) 
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m \ n 2 

In summary, the modified expressions for the virial coefficients are given by Eq. (|l|) for d = 2 and by Eqs. (@) 
and (^7|) (30) for d > 3. We are now in a position to compare the theoretical predictions we have provided in this 
and the previous Section with the available Monte Carlo data for the virial coefficients. We note that the coefficients 
B n yU2 corresponding to n < 5 have been recently evaluated for d — 2 |Lt1JL5| , d — 3 Jf2| , d = 4 jl9j and d = 5 

p9| . Further, the coefficients corresponding to n = 6 for d = 2 are also known jl8|. Because of its special physical 
interest and for reasons that will become apparent later, the case of hard spheres [d = 3) will be addressed separately 
in Section IV. 

We begin our assessment of our estimates of B niy7l2 with the binary mixture of hard discs. This implies using Eq. 
( |32| ) for the original recipe or using Eq. ( |3l| ) in the modified version, which coincides with Whcatley's proposal |33| . 
In Figs. we compare the results of taking these two routes with the ones of Refs. |^0| , p]JT5|Jl^ 1 . As clearly seen in 
the figures, the theoretical prescriptions do a very good job for the whole range of a, the modified version (|l]) being 
almost perfect. 
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FIG. 1. Reduced composition-independent virial coefficients £>* 



m .re — Hi 



from top to bottom, m = 1 and 2. The circles are exact results JlO[, the solid line is the theoretical prediction (|32|) and the 
dashed line is the theoretical prediction 



as functions of the size ratio a for d = 2, n = 3, and, 
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FIG. 2. Reduced composition-independent virial coefficients B£ lin _ ni as functions of the size ratio a for d = 2, n = 4, and, 
from top to bottom, m = 1, 2 and 3. The circles are "exact" Monte Carlo results pH , the solid line is the theoretical prediction 
( [j2~| ) and the dashed line is the theoretical prediction (|3l|). 
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a 

FIG. 4. Reduced composition-independent virial coefficients -B^ ljn _ ni as functions of the size ratio a for d = 2, n — 6, and, 
from top to bottom, n\ = 1, 2, 3, 4 and 5. The circles are "exact" Monte Carlo results ]l8| , the solid line is the theoretical 
prediction (p2j) and the dashed line is the theoretical prediction (pl|). 
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In the case of hard hyperspheres in d = 4 and d = 5, Figs. |5|-[l0| illustrate the performance of Eq. (R) together 
with (|sl)-([il"l) and of Eqs. (0) and (|27|)-(j30|) with respect to the "exact" data of Ref. fl9|| . In these dimensionalities, 
it is clear that while the modified version of the estimates for the virial coefficients continues to exhibit a very good 
performance, the use of Eqs. (^-([TTl) presents limitations for n = 4 and especially for n — 5. Nevertheless, the EOS 
for the binary mixture derived from Eq. (Q) (which yields the original estimates of the coefficients B ni ^ n2 ) 1 once the 
proper Z s (rj) is used in both d = 4 and d — 5, turns out to be very accurate |3q| . 




a 

FIG. 5. Reduced composition-independent virial coefficients B^ i n _ ni as functions of the size ratio a for d = 4, n — 3, and, 
from top to bottom, rt\ — 1 and 2. The circles are "exact" Monte Carlo results [p_9| , the solid line is the theoretical prediction 
(@)~(0) an d the dashed line is the theoretical prediction given by Eqs. (f?|) and (|27|)-([jo|). 
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'0.0 0.2 0.4 0.6 0.8 1.0 

a 

FIG. 6. Reduced composition-independent virial coefficients B* lin _ ni as functions of the size ratio a for d = 4, n — 4, and, 
from top to bottom, ni = 1, 2 and 3. The circles are "exact" Monte Carlo results ljl9|, the solid line is the theoretical prediction 
(@)~(0) an d the- dashed line is the theoretical prediction given by Eqs. (Q) and (p^-(^o|). 
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a 

FIG. 8. Reduced composition-independent virial coefficients B* 1>n _ ni as functions of the size ratio a for d = 5, n — 3, and, 
from top to bottom, n\ = 1 and 2. The circles are "exact" Monte Carlo results [p.9l , the solid line is the theoretical prediction 
(@)~(0) an d the dashed line is the theoretical prediction given by Eqs. (Q) and (p^)-(^o|). 




a 

FIG. 9. Reduced composition-independent virial coefficients B£ ljn _ ni as functions of the size ratio a for d = 5, n — 4, and, 
from top to bottom, n\ = 1, 2 and 3. The circles are "exact" Monte Carlo results illy], the solid line is the theoretical prediction 
(@)~(0) an d the dashed line is the theoretical prediction given by Eqs. (Q) and (p^-(^o|). 
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FIG. 10. Reduced composition-independent virial coefficients n _ ni as functions of the size ratio a for d = 5, n = 5, 
and, from top to bottom, ni = 1, 2, 3 and 4. The circles are "exact" Monte Carlo results [jl9| , the_ solid line is the theoretical 
prediction (fr ) — (|lo[) and the dashed line is the theoretical prediction given by Eqs. (Q) and Q27|)- 



Once we have analyzed the behaviour of the virial coefficients, the important question that now arises is, which is 
the EOS consistent with the modified virial coefficients? In the case d = 2, from Eq. (Blh one gets fl37j 
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(33) 



The EOS consistent with the coefficients (|27|)-(30) can also be explicitly obtained for d > 3 by carrying out the 
resummation of the corresponding virial series. The result is (see Appendix H for details) 
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where K, Kq, K\ and -K2 are independent of the density and are given by 
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K 2 = 



(2 d - 4)(2 2 - d + d - 3)cti: 



{ct? 2 [2 d (d - 2)<7 2 - 8(1 - 2 d ~ 1 + d2 rf - 3 )ai] 
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Here, o\ 2 = (ci + Note that Eq. ( J33| ) for hard discs is included in the structure of Eq. (|34|), except that 

then Eqs. (^5|)-(^8|) cannot be applied and we have instead K = 0, K\ = er 2 (er 2 — (J\), K 2 = o"i(o"i — er 2 ) and 
K = (a 2 ) (1 - XxKxIul - x 2 K 2 /a 2 ). 

Equation j34| ) complies with the exact results given in Eq. (^2|). It is clearly more complicated than the original 
recipe, Eq. ((2^7 and should in principle be useful in particular for mixtures involving components of disparate sizes. 

In contrast to Eq. (||), the proposal ( |34| ) expresses Z m (rj) in terms not only of Z s {rj) but also involves Z s ( and 

We note that if d = 3, it follows from Eqs. (||)-(||) that K = 0, Ki = a 2 {oi - a 2 ) 2 , ^2 = cti(oi - a 2 ) 2 and 
if = (cr 3 ) (l — xii^i/crl — x 2 K 2 /af). Consequently, Eq. (|34|) for hard spheres becomes 



Xl 



1 



1 - 772 



((Tl - <7 2 ) 
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((Tl - (T 2 ) 2 



2(u 2 )u 1 g 2 



It is interesting to note that Eqs. (p3h and (J39|) can be written in the common form 



(39) 



Z m {r\) = Z a (rj) + xi 
1 

+^2 



1 - m 



'/2 



1-% 

- Z B (»?) 



(T 2 — 0\ 



(T 2 
d-1 



d-1 



(40) 



Equation ( ^0[ ) reduces to Eq. ( p3| ) if d = 2 and to Eq. ( p9| ) if d = 3. In addition, it becomes exact for d = 1, i.e., 
Zsiji) — Z m (r)) — (1 — 7 7 )~ 1 - I t se ems th en v ery tempting to assume Eq. ( ptO| ) for arbitrary (i. With such an assumption, 
after using Eqs. JB3[), (B4), (B£) and (B9), we simply get for the composition- independent virial coefficients 



b: 



nior 



n 2 



Tlx 



+ (b ni ~ b n ) —a(a - l)"" 1 + (b n2 - b n ) ^(1 - a) 



\d-l 



(41) 



Of course, this includes the cases d = 2 [cf. Eqj3l])] and d = 3 i c f- E q- ©]• Howeve^Eq. @ differs from the 
prescription (||)-(|^) if d > 4. Although Eq. (§4]) is consistent with conditions (0), @ [for d > 3] and (||), it 
does not reproduce the exact second virial coefficient £?* -Aa) = Vd(l + a) d /2 (except, of course, for d < 3). In fact, 
we have checked that the performance of Eqs. ( |4C| ) and (^|) is very poor for d > 4, so that they will not be further 
considered here. 



IV. THE EQUATION OF STATE OF BINARY HARD-SPHERE MIXTURES 

Up to here, although ocasionally we have mentioned explicit features of hard discs and hard spheres, the develop- 
ments presented in the previous Sections apply for general dimensionality d. Given their intrinsic importance, in this 
Section we focus on the results for the EOS in three dimensions. In order to set the framework of our discussion, we 
shall first recall two other EOS in the literature which share with ours the idea of obtaining Z m (rj) from the knowledge 
of Z s {rj). 

The first one is due to Hamad [^(J who in 1994 |^l[] derived a consistency condition involving the derivatives of the 
contact values of the radial distribution functions gij with respect to J42fl . Imposing this condition, he proposed 
the following equation of state for hard-sphere mixtures 



3*7 

(1-77)3 



T 2\3 



r 3\2 



-(1-77) - 1 



(42) 
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where we have used the label H to refer to Hamad's proposal. It must be noted that Eq. (|4^) was not explicitly written 
in Ref. E^j, where the results are expressed in a rather more involved form. The second EOS for binary mixtures of 
hard spheres that we want to consider was proposed more recently by Barrio and Solana |43] (label BS). It reads 



^ S W = 1 + 7( 1 +^)( 1 + 3 



1 



(a) (a 2 



[ZM - 1] 



(43) 



where (3 is adjusted as to reproduce the exact third virial coefficient ||, namely 

B| = vjlOa 3 , B^ 1 =vl{\ + 2a + ha 2 + fa 3 ) , 



(44) 



More explicitly, in our notation (3 is written as 

= 



B, 



10 



vl(a3)2(l + 3(*)(**)/(*3)) 



(45) 



We now first concentrate on the composition-independent virial coefficients. Our original recipe, Eqs. ([^-([Tol), 
gives for d = 3 
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™2 



(46) 



where we have incorporated the label SYH to distinguish it from the other results. On the other hand, in the 
three-dimensional case, Eqs. (p7|)-(|3"o|) yield 
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(47) 



which coincides with the form proposed by Wheatley p4[ and so the label in this case is W. The virial coefficients 



associated with Eqs. (|42|) and 



are also readily derived. In the former case, first it is easy to obtain 



b n \n{n 1) + 3(n - 1)^^> + |(n - l)(n - 2)g 



2\3 



(48) 



Now, by making m = and m = 1 in Eq. (A2) and m = 1 in Eq. ( |Al|) , and inserting the results into Eq. (^8|), we get 
the composition-independent coefficients arising from Hamad's proposal, 
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(49) 



As far as Eq. (Ea) is concerned, one readily derives (for n > 3) 
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(6„ + /36„_ 1 )( ( 7 3 )"- 1 1 + 3 
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(50) 



Now, after some algebra also involving the substitution of the results obtained by taking m — and m = 1 in Eq. 
(|A2|) into Eq. (p0|), the composition-independent virial coefficients for the Barrio-Solana EOS turn out to be given by 
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(51) 



These four approximations for B* x „ 2 satisfy the requirements ( |17| ) and (|18|). The condition (B2j) is only satished 
by the W approximation, while the condition (^3|) is only satisfied by the H approximation. By construction, the 
BS approximation, E g. (|5l| ) , reproduces the exact third virial coefficients (|44|) . It is interesting to point out the 
approximations (Jhj), ( f47| ) and (|4^) also yield the exact value of the third virial coefficients. 

In order to assess the accuracy of the previous approximate expressions for B* „ 2 with n > 4, we will take the case 
n\ = 3, 712 = 1 as a benchmark for comparison since the coefficient B% x has been obtained exactly for the interval 
< a < 2/V3 - 1 ~ 0.1547 |J. This exact result is 



B 



3,1 



Setting n\ — 3 and ri2 = 1 in Eqs. 
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]), @, @ and (HI), respectively, yields 
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(52) 



(53) 



S*J = ^ [1 + (28 - 6 4 ) a + (5fe 4 - 59) a 2 + 30a 3 ] 



(54) 



and 



B 3 H i = ^[bi-l8 + 9a + 36a 2 + 3 (64 - 9)a 3 ] 



B 
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3,1 



16 L 



15 + (3& 4 - 15)a + (36 4 + 75)a 2 + (96 4 - 45)a 3 ] 



(55) 



(56) 



Note that Eqs. (|53|) and (|5_4|) would coincide if 64 were equal to 18, as happens with the CS equation of state. However, 
they differ if we take the exact value 64 = 18.36477. The maximum deviation of ( |54| ) from the exact result ( p2) is 
about 1.8% (at a ~ 0.08), while that of (||) is about 3.7% (at a ~ 0.09). Also, if 64 were equal to 19, Eq^5|) 

would coincide with Eq. (j54f), while using the exact value for 64 Eq. (|5^) implies that lim Q ^o 

Bf{ = O.Q91u| instead 

which, however, is satisfied by the approximations (p3) and (p4j). Comparison 



3,1 



I,,3 



3' 



of the exact value lim Q _>o B- 

between Hamad's approximation (|5^) and the exact result ([32]) for a < 2/V3 — 1 shows that the relative error of 

(55) monotonically decreases from about 64% at a = to about 18% at a = 2/v / 3— 1. Finally, according to Eq. 

(56) with the exact value of 64, lhn a ^o B^ s — 0.210u|, while in the interval a < 2/\/3 — 1 the relative error of ( |5^ ) 
monotonically decreases from about 16% at a = to about 1.7% at a = 2/V3 — 1. This relative error is always larger 
than that of Eqs. ( f53| ) and but smaller than that of Eq. (j55|). 

Apart from the exact result for B^ l5 as pointed out before, numerical values of the coefficients Z?* n2 corresponding 
to n < 5 are also available for hard spheres |l2|-|l4|,[l6| . In Figs, [ll] and |l2| we display the comparison between the 
theoretical estimates we have just discussed for the virial coefficients and the numerical results. It is clear that the 
estimates Z?*^ „^ , B^ n2 and B^^ 2 for n = 4 and n = 5 are remarkably close to one another and to the exact results, 
except perhaps in the case of B*[ 4 S . On the other hand, the prediction of B^ n is much poorer. 
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FIG. 11. Reduced composition-independent virial coefficients n _ ni as functions of the size ratio a for d — 3, n = 4, 
and, from top to bottom, ni = 1, 2, and 3. The circles are "exact" Monte Carlo results [JD|, the solid line is the theoretical 
prediction (^6|), the dashed line is the theoretical prediction (^), the dash-dot line is the theoretical prediction (^) and the 
dotted line is the theoretical prediction (fsil). 




a 

FIG. 12. Reduced composition-independent virial coefficients -B^ lin _ ni as functions of the size ratio a for d = 3, n — 5, and, 
from top to bottom, m = 1, 2, 3 and 4. The circles are "exact" Monte Carlo results fllrj| , [i"4) , the solid line is the theoretical 
prediction (^), the dashed line is the theoretical prediction (0), the dash-dot line is the theoretical prediction (^) and the 
dotted line is the theoretical prediction (plj). 
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Once we have briefly discussed the performance of the various approximate expressions for the composition- 
independent virial coefficients, we may next wonder how well the approximate EOS for a binary mixture of hard 
spheres perform with respect to the simulation data. To address this issue, we have first to specify what Z s (ij) is. 
Two choices will be made here. On the one hand, since perhaps the most popular and widely known EOS for a single 
component fluid of hard spheres is the CS EOS, we will also use it here. On the other hand, a more accurate EOS 
for the same system is obtained if a rescaled Pade approximation taking into account the known virial coefficients 
is used In Tables [II and IV we present a comparison of the results obtained using the various approximate 

EOS for the mixture with these two choices of Z s (j]), two different mole fractions and a diameter ratio a = 0.3. In 
Table |n|, the values corresponding to the widely used Boubhk-Mansoori-Carnahan-Starling-Leland (BMCSL) EOS 
Pq| have also been included. One immediately notes that, as already stated in the literature |24|,|32|,[55|] , the EOS given 
in Eq. (Q) using the CS EOS for the single component system provides a very accurate account of the simulation data, 
better than any of the other EOS for mixtures considered, including that of Eq. (|39|). If instead of the CS EOS the 
rescaled Pade approximation is taken for Z s (ri), then it is Eq. ( |39| ) which gives in general the best results, although 
the original recipe (||) and the BS EOS ( ff3| ) are not too far behind, the latter being especially accurate in the case 
xi = 0.75. It should be noted that, irrespective of the choice for Z s {rf), both Z^(rf) and Z^ s (r;) for these size ratio 
and packing fractions always underestimate the simulation results, a comment that also applies to the BMCSL EOS. 



TABLE III. Comparison of the compressibility factor from several equations of state with simulation values (Zsimui) for 
binary additive mixtures of hard spheres (d — 3) with a size ratio a — 0.3. The fourth to seventh columns give the deviations 
from Z s i mu \ of the values of the compressibility factor as obtained from the equations of state (^|) , (^) , (0) and ( |43| ) , respectively. 
In these equations of state, the Carnahan-Starling compressibily factor for the one-component fluid |2q] has been used. For the 
sake of comparison, also the results arising from the Boubh'k-Mansoori-Carnahan-Starling-Leland (BMCSL) equation of state 
f|45l are included. 



Xl 
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a 

si mul 


Eq. (2) b 


Eq. 09) c 


Eq. (42) d 


Eq. (43) c 


BMCSL 



0.0625 0.3 2.790 -0.001 -0.009 -0.061 -0.016 -0.014 

0.35 3.473 0.006 -0.011 -0.112 -0.024 -0.020 

0.4 4.410 0.013 -0.019 -0.210 -0.043 -0.035 

0.45 5.722 0.027 -0.033 -0.385 -0.076 -0.063 

0.49 7.158 0.065 -0.036 -0.603 -0.102 -0.081 

0.75 0.3 3.554 -0.005 -0.006 -0.024 -0.010 -0.008 

0.35 4.601 -0.012 -0.014 -0.048 -0.021 -0.018 

0.4 6.045 -0.010 -0.014 -0.079 -0.027 -0.021 

0.45 8.097 -0.002 -0.011 -0.129 -0.033 -0.022 

0.49 10.415 -0.004 -0.021 -0.210 -0.056 -0.037 
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TABLE IV. Comparison of the compressibility factor from several equations of state with simulation values (Z s imui) for 
binary additive mixtures of hard spheres (d — 3) with a size ratio a — 0.3. The fourth to seventh columns give the deviations 
from Zsimui of the values of the compressibility factor as obtained from the equations of state (Q), (^), ( fl2j ) and (|43| ), respectively. 
In these equations of state, a rescaled Pade approximant compressibily factor for the one-component fluid ]24,|44j has been used. 
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V. CONCLUDING REMARKS 



In this paper we have further assessed the merits and limitations of a simple recipe recently introduced f32j| to derive 
the compressibility factor of a multicomponent mixture of d-dimensional hard spheres. We have now considered the 
case of binary mixtures and looked in particular at the composition- independent virial coefficients. The comparison 
with the available exact and simulation results confirmed the usefulness of our approach except perhaps when the 
mixture involves components of very disparate sizes (especially for high dimensionalities). Guided by the conditions 
that the virial coefficients must fulfill in certain limits, a slight modification of the form of the original coefficients 
was made, trying to sacrifice simplicity as less as possible while at the same time improving the accuracy of their 
numerical predictions. It is fair to state that the (modified) composition-independent virial coefficients are in excellent 
agreement with the reported values for all dimensions and all size ratios. 

Due to its particular relevance, the case of the binary mixture of hard spheres (d = 3) was analyzed in some detail, 
both in connection with the virial coefficients and through the consideration of different EOS for this system that exist 
in the literature. As for the virial coefficients, the trends pointed out above about the reliability of the original and 
modified prescriptions hold also for d = 3. T he comparison with the simulation results for the compressibility factor 
of these mixtures given in Tables [II and IV indicates the superiority of the original recipe for the EOS provided one 
uses the CS EOS for Z s (rj). On the other hand, if the more accurate rescaled Pade approximation for Z s (rj) is used, 
it is the new EOS ( j39| ) (i.e., the one obtained after resumming the virial series with the modified coefficients) which 
gives the best performance, although only slightly better than Z m {if) as given by Eq. (|J) with d = 3. Both equations 
are also of comparable accuracy to the one of the presumably best EOS available presently for hard-sphere mixtures, 
namely the rescaled Pade approximant for mixtures introduced by Malijevsky and Veverka |24j. This latter equation 
of state rests on a different philosophy, however, in the sense that it makes use of the known virial coefficients of the 
mixture rather than the Z s (r)), as in our approach and those of Hamad [[lO| and Barrio and Solana E§]. 

In conclusion, we want to add that both the original proposal of the EOS of a binary mixture of d-dimensional 
hard spheres, Eq. (||), and the modified version, Eq. (|34), provide two simple and yet accurate EOS at least with 
respect to the composition-independent virial coefficients and in comparison with the (relatively scarce) simulation 
data. Further confirmation of their usefulness depends on new simulations which we hope will be encouraged by the 
results presented in this paper. 
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APPENDIX A: DERIVATION OF EQS. (@)-(0) 

In order to identify the composition-independent coefficients B ni ri2 from Eq. (^|), we need to expand the moments 
(<r p ). After some algebra, we get 

id^n-Sid-rn^^+Uid-l^din-l) y a d(n- ni -l) (" ~ 3) ! ^ ,- m 

„^ ni!(n-m)! 
x [m(ni - l)(ni - 2)a d + ni(ni - l)(n— n 1 )a d ~ 1 
+ni(m - l)(n - ni )a d - m +n 1 {n- m)(n - m - l)^ - " 1-1 
+ni(ni — l)(fi — ni)a m+1 + n\{n — n\){n — n\ — l)a m 
+ni(n — n\)(n — n± — l)a 

+(n-nx){n-ni-l){n-nx-2)], (Al) 
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+m(n — rii)a m + (n — m)(n — n\ — 1)] . 



where a = 02 /ci- Now we use the identities 
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From here we finally obtain Eqs. (0)-©. 
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APPENDIX B: DERIVATION OF EQS. (34) (38) 



Let us rewrite Eq. (0) as 
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where 012 = (01+02) /2 and the notation 1 <->• 2 indicates that a contribution similar to the three summands appearing 
in the brackets in which the roles of (711,712) and of (01, 02) are interchanged (i.e., ri\ — > n 2 , 01 — > 02 and vice versa) 
should be added. The compressibility factor of the mixture is then given by 
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Now we insert Eqs. (|2TD (|30[) and sum over m and ra. To do that, we use the identities 
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In addition, 
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This allows one to obtain Z m (rj) in terms of Z s evaluated at 77, 7/1/ (1 — 772) and 772/(1 — 771). After some algebra, we 
arrive at Eqs. @)-@. 
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